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Abstract 

On a compact space with non-trivial cycles, for sufficiently small values of the radii of the com- 
pact dimensions, SU(A^) gauge theories coupled with fermions in the fundamental representation 
spontaneously break charge conjugation, time reversal and parity. We show at one loop in per- 
turbation theory that a physical signature for this phenomenon is a non-zero baryonic current 
wrapping around the compact directions. The persistence of this current beyond the perturbative 
regime is checked by lattice simulations. 
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Quantum Chromodynamics (QCD) is the theory of strong interactions. Experimental 
evidence suggests that the theory is invariant under charge conjugation (C), parity (P) and 
time reversal (T) (see pj for a recent account of experimental data). The invariance of 
QCD under P has been rigorously proved in [5]. One of the assumptions of the proof is 
Lorentz invariance, which holds in an infinite volume, but it is manifestly broken at finite 

n 

temperature or in compact space, where parity can be spontaneously broken [3|. Although 
convincing arguments exist ^, a proof of the invariance of QCD under T and C is still 
lacking. 

Recently, it has been pointed out by the authors of that C, P and T are spontaneously 
broken in a geometry with one compact dimension with toroidal topology for sufficiently 
small values of the radius of the torus when periodic boundary conditions are imposed on 
fermion fields. This provides a controllable mechanism for testing the consequences of the 
breaking of those symmetries in QCD. The order parameter is the vacuum expectation value 
(vev) of the Wilson loop winding in the compact direction 

= Tr Pe^^o^^"^^" , (1) 

with a the compact direction of size L, g the coupling and the vector potential. In pure 
gauge S\J{N) {W) oc e'^"" with < n < iV for L < and {W) = for L > L^, where 
Lc is the critical value of the length of the compact direction [6]. Modulo relabeling of 
the axes, the Euclidean rotated system corresponds to the theory at finite temperature and 
the transition that takes places at Lc is the well known confinement-deconfinement phase 
transition [3, Q] . 

When fermions in the fundamental representation are considered, the structure of the ground 
state changes radically. At small radius, if the fermions have antiperiodic boundary con- 
ditions in the compact direction, (W) oc 1 (again this is the case for a system at finite 
temperature in the deconfined phase), while for periodic boundary conditions the Wilson 
loop can take two values with a non-zero imaginary part. These vev are related by complex 
conjugation. Each one of the two values identifies a possible vacuum of the theory. The 
effect of C, P and T is to interchange the vacua. Hence, in this system those symmetries are 



broken. For orientifold gauge theories in the large limit, which are related to QCD jo], on 
a. X S space as the radius of the S is increased above a critical value keeping the radius 
of the small, the system regains invariance under C, P and T 10|. 



The arguments from which the phase structure of QCD on a finite volume is determined are 
based on perturbative calculations. Their validity beyond the perturbative regime has been 



proved by lattice simulations 11 1. 

While the Wilson loop wrapping around the compact direction proves to be useful to charac- 
terise the phases, it is not a quantity that can be accessed directly in experiments. Physically, 
we expect a symmetry breaking to determine a detectable change in the properties of the 
system. Hence, at least one measurable quantity that is not invariant under the broken 
symmetries should acquire a vev. The spatial components in the compact directions of the 
baryonic current ji = X]n=^i "^"^ ® lii^n, where ipn is the fermion field for flavour n, the sum 
runs over the flavour index and I is the identity in colour space, satisfy this requirement jl^ . 
Moreover, like the system, they are invariant under CP, CT and PT. This makes the ji suit- 
able candidates as detectors of the symmetry breaking. If (j^) 7^ for the compact direction 
a, an observer will see a non-zero flux of baryons in that direction. 

Using a similar ansatz to the one of [5!|, we shall now show at one loop in perturbation theory 
that indeed the vev of the spatial current in a compact direction is different from zero. The 
Lagrangian for a SU(A^) gauge theory coupled with Nf degenerate flavours of fermions of 
mass m in the fundamental representation is 

Nf 

C = --^Tr {G,,{x)G^^''{x)) + M - 4 - M^) , (2) 

n=l 

where G^i, = d^A^ — d^A^ + [A^, Ay]. The corresponding partition function is 
Z = j (VA,) det ^ - 4 - mf^- ^ ^'-^no.^<^^n 

We consider the system on a T x L'^ manifold, in which the T direction corresponds to 
time and the three spatial compact directions L are equal. We impose periodic boundary 
conditions in space, while T is assumed to be large enough for the choice of boundary 
conditions in that direction to be irrelevant. The path integral ([3]) can be evaluated at one 
loop, by fixing a diagonal background gauge. This gives an effective one loop potential for 
the diagonal components of the gauge field 



A 



L 



^Vi = Q mod(27r) , (4) 



which reads 12 



V{vi, ...,vn) 



■ N N 

J2fiO,v,-v,)-2NfJ2fi 

.i,j = ^ 



m,Vi) 



i=l 



(5) 



The first sum comes from the integration of the fluctuations of the gauge and ghost fields, 
while the second sum comes from the fermion determinant. The function / is defined as 



K2{mLk) . 2 Ar ^ 
I , 7T, sm \ -k ■ V 



(6) 



with the sum running over vectors in Z'^ — and K2 the order two modified Bessel function 
of the second kind. From the asymptotic behaviour of K2{x) at small x, K2{x) = 2/x'^, we 
get 



f{0,v) 



E 

fc^O 



sm 



2/1 



(7) 



For m > K2{mLk) ^ e"^^^ tx / {2mLk) and the sum in / is dominated by terms 
with k = 1. This is true also in general, the higher frequencies in the sum being quickly 
oscillating with amplitude suppressed at least as l/k"^. With the constraints in Eq. (jl]), the 
minima of the effective potential are located at 
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■ N 



TT 



TT 



for odd 
for N even 



(8) 



There are eight degenerate minima for odd and one minimum for even N. In the former 
case, (W) develops a vev with an imaginary part, and the spontaneous symmetry breaking 
occurs. 

The baryonic current can be computed adding a source to the Lagrangian ([2]). Defining 



C{fl) = C + fl-j , 



(9) 



we obtain 



(10) 



fl=0 
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with Z[p] the partition function in the presence of a source /2. 

The source has the effect of shifting Vi^Vi + Lfl in the expression for the effective poten- 
tial ([5]). This does not change the gauge contribution. Since 

at the minima (El) we get 

NfN fmLY ^ K2(mLk) . fr A r 

MO 

(ji) is zero when Vi = or Vi = n (i.e. when the symmetry breaking does not occur in 
direction i), is odd under Vi — fj, goes to zero when m —>■ oo. Hence it fulfills all the 
natural requirements in the current scenario. In particular, we expect a non zero current for 
an odd number of colours. 

In order to get a better handle on the properties of the baryonic current in the broken phase 
beyond perturbation theory, we have performed a lattice simulation using four flavours 
of staggered quarks coupled to an SU(3) gauge field. The number of flavours has been 
flxed as the minimal one for which the staggered action has an undoubtedly well deflned 
continuum limit. For the pure gauge action we have used the standard Wilson form Sq = 
(3 (l — |Trf/p), where /? = IN/ is the coupling of the theory. Up is the path-ordered 
product of link variables around the elementary plaquette P and the sum runs over all 
plaquettes P. For the fermionic part we have used the simple staggered action 

'S'f = J^?7^(a;)^ (x(a;)f^A.(2;)x(a; + A) - c.c.) + am J^x(a;)x(a;) , (13) 

X,^ X 

with ?7^(x) = (— l)Et;=o {rio{x) = 1), x a complex three- vector, am the mass in lattice units 
(a is the lattice spacing) and c.c. stands for the complex conjugate term to the flrst one in 
parentheses. More complicated formulations of the action or choice of another discretised 
form for the fermionic fields would have added extra complication with very little payback 
for the problem at hand. 
Using as a base the publicly available MILC code we have performed a simulation for 
/3 = 5.5 and am = 0.1. The physical scale has been determined by measuring the Sommer 
parameter tq jl^ on a 24 x 16^ lattice, where the three equal spatial directions Ng have been 
closed with periodic boundary conditions and the temporal direction Nt with antiperiodic 
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FIG. 1: Scatter plot for 1000 measurements of the Wilson line in one compact direction on a 24x4^ 
lattice at P = 5.5 and am = 0.1. The directions corresponding to the three cubic roots of the unity 
are indicated by the black solid lines. 

boundary conditions for the fermions, while the gauge fields are periodic in all directions. 
We find aro = 4.0(1); since the Sommer scale is ~ 0.5 fm, the lattice spacing is a ~ 0.125 
fm, which means that = aNg ~ 2 fm and Lf = aNt — 3 fm. Hence, in physical units the 
lattice is large enough for the calculation to be reliable and the spatial volume is such that 
C, P and T are not broken. 

We then studied the system with the same (3 and m on a 24 x 4^ lattice, with the same 
boundary conditions as above. The spatial geometry is a three-torus, while the size of the 
temporal direction is large enough for the system to be confined. In this setup, Ls — 0.5 fm. 
A quick check of the vev of the spatial Wilson loops shows that the system is in the broken 
symmetry phase. An example of the obtained distribution for (W) is displayed in Fig. [T], 
which shows (W) clustering around e*3'^. 

The discretised version of the baryonic current can be obtained like in the continuous case, 
by relating the physical fermionic degrees of freedom to the staggered ones. Defining the 
massless Dirac operator as 

^"'^ = E^A<(^) (f^M(^) W - Uli^ - i^)Sy,.-f^) (14) 
and the four matrices 

= {U^ix)5y,^+f, + Ulix - (15) 
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FIG. 2: The imaginary part of the current and the imaginary part of the Polyakov loop in one 
compact direction as a function of the Monte Carlo sweeps. 
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FIG. 3: As in Fig. [21 but in another compact direction. The system shows a transition between 
two vacua. The transition probability is finite, due to the finite lattice extension. 

the current reads 

(j) = ^Tr((D + m)-^i?) . (16) 

In order to evaluate the current on a given configuration, we have taken 100 stochastic 
estimates. Since the current is an antihermitian operator in the Euclidean space, we expect 
its imaginary part to develop a vev, while the real part should average to zero. In Fig. [2] 
we show the behaviour of the imaginary part of the baryonic current in a compact direction 
as a function of the Monte Carlo sweeps, and we contrast such behaviour with that of the 
imaginary part of the Wilson line in the same direction. Not only does the plot show that 
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the baryonic current is different from zero, but it also strongly suggests that there is a 
correlation between the value of the current and the value of the Wilson line. In particular, 
the modulus of the imaginary part of the current grows when the modulus of the imaginary 
part of the Wilson line grows, the sign being opposite between the two. This is better shown 
by Fig. [3l which displays the behaviour of the current in another compact direction. In 
this case, the system makes a transition between the vacuum identified by the phase of the 
Wilson line being ^tt to the other vacuum and then back. Noticeably, the current changes 
sign exactly at the points in which the imaginary part of the Wilson line changes sign, with 
its magnitude always tracking closely the magnitude of the phase of (W). The sum of the 
terms with = 1 in the current (fT2i) is proportional to (W). The strong correlation between 
the two quantities suggests that the non-leading terms in f|T2|) do not affect significantly the 
behaviour of ja- 

Transitions between the different vacua like those shown in Fig. [3] are possible because 
of the finite spatial size of the system. We have verified that increasing (3 at fixed lattice 
extensions A^^ and Nt, which corresponds to decreasing the physical volume, the frequency 
of the transitions increases. Likewise, decreasing j3 decreases the likelihood of a transition 
taking place. 

Since the baryonic current is zero for symmetry reasons in the symmetric phase, its behaviour 
in the broken symmetry phase makes it legitimate to use that current as an order parameter 
for the symmetry breaking. For consistency, we have also checked that the real part of 
the current in the compact directions and the zero component are zero also in the broken 
symmetry phase. 

In order to evaluate the magnitude of the current, we averaged over directions for which no 
tunneling between the two vacua took place. We find 

|Im I = 0.060 ±0.002 . (17) 

It is instructive to compare this number with the one loop expression, Eq. (fT2l) . which gives 
(jq) ~ 0.037473(4), where the error is a conservative estimate for the truncation of the sum. 
Hence, quite remarkably the one loop calculation pins down the correct order of magnitude 
even for a compact dimension with size of the order of 1/Aqcd- Non-perturbative effects 
could explain the discrepancy between the perturbative formula and the measured value. 
Besides, our calculation being at one single lattice spacing, we do not have any handle 
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on the size of discretisation errors. For this reason, a careful comparison between the 
perturbative expression and the lattice result should be the subject of a more detailed 
study, which is beyond the scope of this paper. Our preliminary Monte Carlo results for the 
current closer to the continuum limit show substantial agreement between the measured 
value and the perturbative formula. 

In conclusion, we have shown that QCD on small compact dimensions with non-trivial 
cycles is characterised by a flow of current whose sign depends on the vacuum selected by 
the system. The persistent baryonic current reminds the supcrcurrent observed in super- 
conductors. However, there is a fundamental difference: unlike the case of superconductors, 
in QCD in compact not simply connected space the current is still conserved, since the 
U(l) baryon symmetry (which in the case of QCD is a global symmetry) remains unbroken. 
The persistent flow is induced by the spontaneous breaking of a discrete symmetry, charge 
conjugation. The baryonic current can be used as an order parameter for the spontaneous 
breaking of charge conjugation in SU(A^) gauge theories. Over the Wilson line, it has the 
advantage of being an observable quantity. Moreover, unlike the Wilson line, which is 
ultraviolet divergent on the lattice, the baryonic current is a well defined observable. This 
makes it better suited for numerical studies of the physics of C parity spontaneous breaking 
close to the continuum limit. A similar investigation is currently in progress. 
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